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We attempt to explain the large A/p ratios measured in heavy ion collisions at 12 A • GeV/c 
beam momentum within an hadronic framework. The ratios are large compared to corresponding 
ratios in pp collisions, and to thermal fits. We show using a simple model and a detailed cascade 
calculation that different annihilation cross-sections of A's and p's, and the net conversion of p's to 
A's, do not account for the enhancement. Uncertainties in elementary cross-sections and formation 
times are also considered. 



I. INTRODUCTION 



Anti-Lambda (A) cross-sections in heavy ion collisions are of interest because strangeness production is a potential 
signal for QGP formation. The ratio A/p is of special interest since it reflects the production of s-quarks 

relative to non-strange light anti-quarks, and should increase substantially relative to the production expected from 
a superposition of NN collisions, if a QGP is formed. 

Recent experiments have reported [^-Q measurements of p and A production in various heavy ion systems at the 
Brookhaven AGS and at the CERN SPS. Experiment E859 reports the ratio of A and p rapidity distributions to 
be 3 ± 1 ± 1 in central Si+Pb collisions. Q This ratio is corrected for finite experiment acceptance, efficiencies, and 
for p creation from A decay (feed-dovirn). It also takes into account that neutral S particles cannot be distinguished 
from the A sample. Experiment E864 has measured § the p production cross-section in Au+Pb at 11.6 A • GeV/c. 
They compare this measurement with a similar one from the E878 collaboration that was obtained with a focusing 
spectrometer, and interpret the difference between the two measurements as an indicator of A production. E864 
estimates that A/p > 2.3 at the 98%CI, with a most probable value of 3.5. At the SPS, NA35 has pubhshed [|j7) A/p 
ratios for pp, pA, S+S, S+Ag, and S+Au collisions at 200 A • GeV/c, and they observe a significant rise from 0.25 for 
pp-collisions to 1.5 for the heavy ion systems. 

It is tempting to interpret the large reported ratios as evidence for the formation of a QGP: Many authors ]^ [ic|] 
have argued that the CERN multi-strange baryon ratios pl|-p^ can only be described by a QGP scenario. This 
conclusion is, however, challenged by E864 Au+Au RQMD simulations that find that the p production cross-section 
appears to be lower than expected from the scaled NN p cross-section. 

Here, we shall try to address quantitatively all possible hadronic contributions to the A/p ratio. We restrict 
calculations to AGS energies and show that "differential annihilation" of the two species and p-to-A conversion 
processes can indeed enhance the A/p ratio, but we conclude that the effect is not large enough - thus hinting at a 
production mechanism outside of the standard hadronic interactions. We present these arguments as follows: The 
observed A/p ratios cannot be explained by a thermal model, unless severe inconsistencies with other measured data, 
such as the charged pion and kaon cross-sections, are introduced. This has been shown by many authors, [l^-|l^] 
and we give our own calculation in the next Section. Thermal equilibration is unlikely at present energies, so we 
consider transport simulations in the rest of the paper, first discussing the relevant cross-sections, then the results 
for a simple geometrical model, and finally the results of a detailed cascade calculation. Actual measurements of 
the A and p cross-sections [p^p^ are used, not event generator parametrizations. We shall consider uncertainties 
in these calculations, particularly those arising from the relatively poorly known A annihilation, and from particle 
formation times. Speculative conclusions based on the quantitative discrepancy between calculations and the actual 
measurements are presented in the last Section. 



*Submitted to Proceedings of the 14th Winter Workshop on Nuclear Dynamics, Snowbird, Utah, 31 January - 6 February 
1998, "Advances in Nuclear Dynamics 4," W.Bauer and H.G. Ritter, eds. (Plenum Publishing, 1998). 
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II. THERMAL MODEL 



In a thermal model, the final relative abundance of a particle is determined by both the primary number of this 
species at freeze-out, and by feed-down from heavier species after freeze-out. We report here results for central 
Au+Au collisions at 11.6 GeV/c, and include all mesons with rest mass < 1 GeV/c^, and all baryons with rest mass 
< 1.7 GeV/c2. 

We assume thermal and chemical equilibration for freeze-out, so that all relative abundances can be obtained from 
four parameters: the freeze-out temperature Tq, the electric chemical potential fie, the baryonic chemical potential 
/ib, and the strange chemical potential Hs- Here, 



Si fJ,s 



(2.1) 



are particle chemical potentials, where qi, hi and Si are the charge, baryon number and strangeness of species i, 
respectively. The thermal fit parameters are obtained by applying the conditions \QI B — 0.40| < 1% and IS*] < 2%, 
together with the constraints listed in Table 1. We purposefully choose rather large ranges in these latter constraints 
to show how difhcult it would be to obtain the experimental K/p ratio.^ 

In Table 1, the "errors" on the best fit values for Tq and the chemical potentials are ranges that lead to results 

consistent with the constraint intervals.^ We see that A/p^l.9, at best. It is instructive to consider the results 
graphically. Figure 1(a) shows K/p as a function of K'^ /K~ for various freeze-out temperatures Tq. Large K/p > 2 
result only if the freeze-out temperature and/or observed / ratio are pushed unreasonably high. A similar 
conclusion follows from Figure 1(b), which shows A/p as a function of if "•"/tt"'", with Tq = 120 ± 14 MeV. 

We conclude that the experimental ratio K/p ratio has, at the very least, a significant non-thermal component: the 
lower experimental bound for the ratio reported by E864 is larger than any reasonable thermal fit would allow. 



TABLE 1 Thermal vs experimental particle ratios for central Au+Au collisions at 11.6 GeV/c. The parameter 
ranges are To = 120 ± 14 MeV, fib = 556 ± 19 MeV, fi, = 111 ± 14 MeV, and fie = -14 ± 2 MeV. 



Ratio 


Constraint 


Thermal Model 


Data 


Ratio 


Rapidity Ref. 


K+/7r+ 


0.16-0.28 


0.23 ±0.03 


0.22 ±0.01 


0.5-1.3 [ 


L9 




K+/K- 


4.0-6.0 


4.73 ±0.53 


5.0 ±1.0 


0.5-1.3 [ 


12 




K-/TV- 




(3.50 ±0.62) X 10"^ 


0.028 


1.2-2.0 1 


20 




n+/p 


0.6-1.2 


0.71 ±0.09 










0.8-1.4 


1.00 ±0.10 


1.00 


1.2-2.0 1 




A/p 




0.16 ±0.02 








p/p 




(3.48 ± 3.44) X 10"* 








A/p 




1.58 ±0.30 









III. CROSS-SECTIONS 



Generally, a thermal description is useful for a given particle species if its mean free path is small compared to 
the system size. This fact alone means that other approaches should also be investigated. We do so in the next two 
sections, and discuss here the most important physical input, viz., the relevant cross-sections. The following are the 
relevant processes we shall consider: Firstly, the production of p's and A's: 

bi+b2^ N1 + N2 + B + B i+n) (3.1) 
M + b^N + B + B (±7r) (3.2) 
M + M ^ B + B (±7r) (3.3) 
S+ + S- ^ B + B (±7r) (3.4) 

(3.5) 



^Our number of A's throughout includes E^'s, as experiment does not distinguish between the two species. 
^We note that these parameters are not inconsistent with overall energy conservation considerations. Recall also that mt 
spectra are blue shifted by the presence of strong flow. 
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(3.6) 
(3.7) 



Here, b represents a non-strange baryon, N a nucleon, B any baryon, M a light unflavored meson, and a ±1- 
strangeness meson. Secondly, (same notation): 



p + b 
b + A 

S+ + b- 

M + 

M + A- 



X 

A + TT 
A + TT + 

b + S+i+7:) 



(3.8) 
(3.9) 
(3.10) 
(3.11) 
(3.12) 



We refer to ( |3.8|) and (3^) as annihilation processes, and ( 3.10 )-( 3.12 ) as conversion processes. 

The experimental pp annihilation cross-section is well known, and can be para- 
metrized as 



\Plab) 



67 pj' mb 



(3.13) 



where piab is the momentmn of the "beam" particle in GeV/c with the "target" at rest. The solid line in Figm'e 2 shows 
this parameterization of the data (diamonds). The A annihilation cross-section, on the other hand, is relatively poorly 
known, especially in the energy range we are interested in. We model it by assuming that the elastic cross-sections 
for pA and pA are equal, and then use the data of Ref. |lq| to obtain: 



anm, 



^{Plab) = 15 (p,a6/10)-" mb , 



(3.14) 



where the same comments apply as for Eq. ( p. 13^). While the data is best fit by a = 0.5, the uncertainty is rather 
large (see Figure 2; the stars are data from Ref. |18|). In fact, one might well argue that the data is consistent with 




0.15 



FIG. 1. (a) Thermal A/p ratios as a function of /K~ for various freeze-out temperatures To. (b) Thermal A/p ratios as 
a function of A'^/tt^, with To — 120 ± 14 MeV. Error bars and scatter points indicate values consistent with constraints not 
shown; the lines are to guide the eye only. 



3 



400 




I I \ , , L 

4.00 lO.OO 40.00 

Pia.b [GeV/c] 



FIG. 2. The annihilation cross-sections of p (A) with nucleons, as a function of the incident momentum of the p (A). 

a = 0.^ For a = 0.7 the A data is practically indistinguishable from the p data. We shall subsequently investigate 
the behavior of A/p with a. 

The p and A scattering with mesons is also an important process we need to consider. Broadly speaking, we have 
three types of collision: (1) thermalization of p's and A's through elastic collisions; (2) production of resonances that 
eventually decay back into p's or A's; and (3), most importantly, net conversion of p's to A's. Chief amongst these is 



P 



+ — > A (or resonances of A) 



(3.15) 



for which we know the charge conjugate reaction to have a sizeable cross-section. The process ( 3.15| ) thus contributes 
significantly to reducing the p abundance while enhancing the A abunda nce i n the final state. Given the pronounced 
strangeness enhancement in large systems such as Au+Au, the process (3.15) should be a relatively important piece 
of the A/p "puzzle." 

Finally, we shall need the A/p ratio in pp colfisions. At 20 GeV, it has a value of 0.25-0.30. At AGS 

energies {y/s ~ 5 GeV) its value is less established. Using Refs. pl|-p3[, we infer a value of ^ 0.2, but we shall use a 
value of 0.25 throughout this work. This implies that the heavy ion A/p ratios we obtain in subsequent calculations 
will be upper bounds. 



IV. GEOMETRIC MODEL CALCULATION 

It is often rather useful to have a simple understanding of cascade output. We shall perform such a detailed cascade 
simulation in the next Section; here we model p and A production, annihilation and net conversion {p, n to A, S) in 
a simple geometric model. Our starting point are the rate equations 



dNr, 



1 



1 



Np dz 



1 



Aa 



N^dz + — Np dz 



(4.1) 



where Aa and Ac are the annihilation and net conversion mean free paths for p's, respectively, and A^ is the annihilation 
mean free path for the A. We take the initial A/p ratio from pp collisions: 



^QCD sum rules imply a behavior consistent with a ~ 0. 
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0.25 



as discussed in the previous section. The final number of p's or A's depends on the distribution of the length z of 
nuclear matter that the anti-particle passes through: 

N = {N{z)) . 



Integrating Eqs. (4.1) gives the number of p and A after they have passed through a fixed length z of nuclear matter: 



Np{z) = NP e-(i/^»+i/^=)^ (4.2) 

Nj,{z) = < e-^/^" + 2 ("e-^/A. _ ^-{i/x.+i/x^)A _ (4^3) 

-'^ AT ^ AT ^ / 

The factor of two has been introduced to account for the net conversion of n's to A's. 

Next, we need to model the geometry of the collision. We shall assume that z is related to the combined thickness 
of the beam and target nuclei t via 

z ^ I3t , 

where /3 < 1 is a constant. It measures the effect of matter expansion, the local average momentum distribution, etc. 
The distribution of t in a central AA collision is 

m - ^ , (4^4) 

where R is the nuclear radius. This equation allows us to compute the survival probability Ps = (e~^/^) of a p or A 
as a function of R, A and /3. The analysis for AB collisions is similar. 
It remains to discuss the mean free paths. Generally, 



T- = ( dpiab PiiVlab) CTji{piab) ^^r~^ ' 

Aj ^ J dpiab 



(4.5) 



where piab is the momentum of the p or A measured in the rest frame of the target i (nucleons, pions or kaons). Also, 
Pi is the density of target species i, ajt is the corresponding cross-section, and dNj/dpiab is the distribution of piab- 
Since the factor (3 already takes into account any changes in density, pi should be regarded as the density for which 
no expansion occurs. 

A rough approximation for Aa follows assuming that the p and A are produced at rest in the NN frame, and are 

annihilated by the nucleons streaming by. In this case, piab ~ \J ^ '^■^ GeV, with p ^ 0.16 fm""^, and 

therefore Aa ~ 1.6 fm.^ A more detailed calculation that accounts for residual p motion in the NN frame gives 
Aq ~ 2.0 fm. Similarly, using appropriate cross-sections and thermal densities, we obtain 

A„ w 4.2 (3.76)"" fm 
A(i) « 180 fm 
A(2) w 88 fm 

where A^^'' and A^^^ are the net conversion mean free paths due to collisions with (N,7r) and K~^, respectively. 

The dependence of the two components of the A/p ratio {viz., annihilation and conversion), are shown in Figure 3 
for a central, symmetric {AA) collision, as a function of the effective size of the system, i?e/ / = PR. The ratio 



*Note that z/\ is Lorentz invariant; we perform calculations in the rest frame of the nucleus. 
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FIG. 3. A/p as a function of / = I3R for a central AA collision, with a = and zero formation time (solid line) . Individual 
contributions from annihilation and conversion are also shown. 

saturates for large i?e//^6 fm/c because peripheral collisions play a increasingly important role; the ratio of 
survival probabilities tends to the ratio of mean free paths cubed. On the other hand, for truly infinite matter all 
anti-particles pass through a common large length, and A/p is exponentially increasing with Reff- 

The significant discrepancy between the data and our calculation (see also next Section) may indicate that the 
ratio of core-produced p's to peripherally produced p's is larger than predicted by the simple binary hadron-hadron 
collision scenario. 

The solid lines in Figure 4 show the final values of K/p from this simple geometric model, if we use (3 = 0.5. The 
agreement with the cascade calculation (circles; see next Section) is very good. 




V. CASCADE CALCULATION 



The production rate of p's or A's is only O(10~^) per event at the AGS, and cascade calculations are rather CPU 
intensive. We shall therefore perform an effective calculation of the survival probability by putting p's or A's in by 
hand, wherever and whenever a collision occurs in which the energy is sufficient for a pp or AA pair to be produced. 
Thus we assume that the pair production does not depend on energy, once above threshold.^ The evolution of 
nucleons, pions and kaons is not allowed to be influenced by the presence of p's or A's - we restore particles that 
interacted with p's or A's to their pre-collision kinematics. 

We may then calculate the A/p ratio via the survival probabilities Pg 

^^.,f^^ + 2.4, (5.) 

P Ps[p) p" 

where p" is the number of surviving p's and A'^ the number of converted A's. Also, 0.25 ± 0.1 is the A/p ratio in 
pp collisions, while Vc ~ 0.25 is a correction factor.^ The factor two accounts for the conversion of n's to A's. 

To illustrate the effect of differential annihilation versus net conversion, we show in Table 2 the different components 
of Eq. ( p3.lD for different a and r = in central Au+Au collisions at AGS energies. For small a (small A annihilation 
cross-section), A/p enhancement results almost entirely from differential annihilation. As a increases, the survival 
probability of a A and a p become equal, and the sole enhancement in the final ratio results from conversion. 



^Also, we assume that the production is not medium dependent. 

^We overestimate A-production because we assume that all collisions of non-strange anti-baryons with positive strange mesons 
result in a A. 
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TABLE 2 Survival probabilities, the conversion rate, and final A/p 
ratio for various a and r = fni/c, in a central Au+Au collision. 



a 


P,(A) (%) (%) AVp= 


A/p 


0.0 
0.2 
0.4 
0.7 


12.1 ±0.3 1.5 ±0.1 0.75 ±0.07 
6.0 ±0.2 1.5 ±0.1 0.45 ±0.06 
3.3 ±0.2 1.5 ±0.1 0.38 ±0.06 
1.5 ±0.05 1.5 ±0.1 0.16 ±0.02 


2.4 ± 1.0 
1.2 ±0.5 
0.7 ±0.3 
0.3 ±0.1 




FIG. 4. Cascade A/p ratios for various a, t and systems. The solid lines are for the geometric model with /? = 0.5. 

We show in Figure 4 the final A/p ratios for several values of a and t for various systems. The solid lines are results 
from our simple geometric model using /3 = 0.5, and show that the simple geometrical model predicts the ratios rather 
well. We note that the ratio is larger in Si+Pb collisions than in Au+Au collisions, both in our calculation and in 
experiment. Also, note that (1) realistic values for a are probably closer to 0.5 than 0.0; (2) formation times are 
probably not zero; and (3) our initial {i.e., pp) A/p ratio is probably an over-estimate. We thus conclude that in our 

present calculation A/p ^ 1 and most definitely ~ 2 for the Au+Au system (~ 3 for Si+Pb). 



VI. CONCLUSIONS 

The A/p ratio in NN collisions at AGS energies is ~0.25, heavy ion experiments measure much larger values. 
The ratio can only be pushed to values above 2 with much difficulty in a hadronic thermal model with chemical 
equilibration. For Au+Au, using upper limits of /-k^ = 0.30, / K~ = 6.0, and Tq = 140 MeV, we obtain an 
upper limit of A/p ^ 2, below the E864 98%-lower confidence limit of 2.3. 

We have considered a non equilibrium description, and presented results from a simple geometric model and a 
detailed cascade calculation. The geometric model does rather well in reproducing trends in the data. Chief inputs to 
our calculations are the A and p annihilation cross-sections with nucleons, as well as processes that convert p's into 
A's. The A annihilation cross-section is relatively poorly known, but it may be argued that it is somewhat less than 
p annihilation at AGS energies. Our parameter a controls this difference, with a = corresponding to the greatest 
difference, and a = 0.7 corresponding to a cross-section that is practically indistinguishable from p annihilation. 
Currently, the cross-section data is best fit by a = 0.5. 

In Au+Au cascade simulations at AGS energies, the largest A/p ratio obtainable is ~ 3.5, for a = and r = fm/c. 
Both these input values are extreme. What we currently believe to be more reasonable inputs lead to a value of 
A/p ~ 1, far from the lower bound of the E864 experiment. 
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Interestingly, higher A/p ratios result in central Si+Pb collisions, both in our calculations and in experiment. Within 
our model this has a simple geometric explanation: For Si+Pb, the p's and A's are produced close to the beam axis, 
leading to a larger average nuclear thickness that the anti-particle must traverse. We also note that at higher energies 
the difference between the A and p annihilation cross-section practically vanishes. One might therefore predict (in 
hindsight) that A/p at the SPS should be lower than at the AGS, as indeed observed. 

We conclude that the large A/p ratios in AGS heavy ion collisions are not easily explained by hadronic mechanisms. 
QGP formation might be a possible solution, but a more accurate measurement of the A annihilation cross-section in 
the relevant energy range is sorely needed before any more definite conclusion can be reached. 

This work was supported in part by the U.S. Department of Energy under Grant No. DE-FG02-93ER40713. 
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